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ABSTRACT 

This paper presents an analytical simulation of the dynamics of a gear 
transmission system compared to experimental results from a gear noise test 
rig at the NASA Lewis Research Center* The analytical procedure developed 
couples the dynamic behaviors of the rotor-bearing-gear system with the 
response of the gearbox structure. The modal synthesis method is used in 
solving the overall dynamics of the system. Locally each rotor-gear stage is 
modeled as an individual rotor-bearing system using the matrix transfer tech- 
nique. The dynamics of each individual rotor are coupled with other rotor 
stages through the nonlinear gear mesh forces and with the gearbox structure 
through bearing support systems. The modal characteristics of the gearbox 
structure are evaluated using the finite element procedure. A variable time 
stepping integration scheme is used to calculate the overall time transient 
behavior of the system in modal coordinates. The global dynamic behavior of 
the system is expressed in a generalized coordinate system. Transient and 
steady state vibrations of the gearbox system are presented in the time and 
frequency domains* The vibration character istics of a simple single mesh gear 
noise test rig is modeled. The numerical simulations are compared to experi- 
mental data measured under typical operating conditions* The comparison of 
system natural frequencies, peak vibration amplitudes, and gear mesh frequen- 
cies are generally in good agreement. 
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NOMENCLATURE 

modal function of the i th mode in x-direction 

modal function of the i th mode in 8-direction 

modal function of the i th mode in y-direction 

gearbox damping matrices 
torsional damping matrix 

bearing direct and cross-coupling damping matrices 
bearing excitation forces 
gear mesh torque 

gear mesh force in x- and y-directions 
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external excitation moment 
external excitation forces 
gryoscopic-angular acceleration matrix 
gyroscopic-angular rotation matrix 
identity matrix 

rotational mass moment of inertia matrix 
average stiffness matrix 
gearbox stiffness matrix 

compensation matrices in x- and y-direction 
shaft stiffness matrix 
torsional stiffness matrix 

gear mesh stiffness between i th and k th rotor 
bearing direct and cross-coupling stiffness matrix 
mass-inertia matrix of rotor 
mass-inertia matrix of gearbox 

the number of modes used to define each motion 
radius of gear in the i th rotor 
gear generated torque 

generalized motion in x- and y-directions 
gearbox motion in x-, y- and z-directions 
gearbox motion at bearing supports 

gear displacements in x- and y-directions of the i th 
rotor 

gear forces in x- and y-directions 

motion of rotor at bearing support 

angle of tooth mesh between k th and i th rotor 

rotational displacement of the i th rotor at the gear 

rotational displacement vector of the i th rotor 
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lateral and torsional eigenvalue diagonal matrices 


[A 2 ],[A t 2 ] 

[$] , [$ ] lateral and torsional orthonormal eigenvector matrices 

k t k 4-h 

of the k rotor 

<f> j th orthonormal mode of k th rotor at 1 th station 


INTRODUCTION 

The study of gearbox vibration is an important area of research among the 
engineering community* With increased speed and torque requirements and 
reduced allowable vibration levels, a critical need exists to develop an 
accurate analytical simulation of the dynamics of gear transmission systems. 
Recently, a number of analytical modeling and simulation procedures have been 
developed to predict the dynamics of multistage gear transmission systems. 

Some experimental work has also been performed to improve the understanding of 
gear stress and vibration. There exits a need to correlate analytical and 
experimental efforts so that the analytical methods can be refined and 
verified . 

There is a wealth of literature on vibration analysis and simulation of 
gear transmission systems. August (1986) simulated the vibration characteris- 
tics of a three-planet transmission system. Pike (1987) and later Choy 
(1988b) used the gear relationships introduced by Cornell (1981) to study the 
dynamics of coupled gear systems. Mitchell (1987) applied and matrix transfer 
method to simulate gear vibrations. Ozguven and Houser (1988) and Kahraman 
(1990) used a finite element model to predict the dynamics of a multigear mesh 
system. Choy (1991) applied a modal synthesis technique in conjunction with 
finite element and the matrix transfer methods in both time and frequency 
domains to calculate the transient rotor and casing motions. 

There are few studies correlating analytical predictions and experimental 
results. Some experimental correlation work by Lim (1990) compared analytical 
predictions of vibration of a gear housing and mounting with measurements. 

This paper reviews the development of a global dynamic model of a gear 
train system, and presents results of a comparison made between the analytical 
predictions of the model and experimental results from a gear test rig. A 
combined approach using the modal synthesis and finite element method for ana- 
lyzing the dynamics of gear systems was developed. The method includes 
dynamic coupling between the housing and the gear-rotor system. The single 
mesh gear noise test rig at NASA Lewis was used as an example in this analy- 
sis. The gearbox modal characteristics were evaluated by solving a finite 
element model. These modal characteristics were compared to those measured in 
the experimental study to insure accurate modeling of the experimental rotor- 
gear-housing system. Transient vibrations of the rotor-gear stages and the 
gearbox housing were analytically evaluated by the modal method. A numerical 
FFT (Fast Fourier Transform) algorithm was used to examine vibration in the 
frequency domain. Frequency spectra of the predicted gearbox vibration were 
compared to those measured in the experimental study. The following discus- 
sion and conclusions are drawn from these comparisons between the experimental 
results and analytical predictions. 
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ANALYTICAL PROCEDURE 


The equations of motion for a single mesh multimass rotor-bearing-gear 
shaft system, with the coupling effects of gearbox vibrations and the rotor 
inertia-gyroscopic effects, can be written in matrix form for the i rotor 
(Choy 1987, 1989) for the X-Z plane as: 
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Here F and F are force excitations from the effects of mass imbalance 
and shaft residual bow in the X- and Y-directions . and F are the X 

and Y gear mesh forces induced from gear teeth interactions. The bearing 
forces are evaluated from the relative motion between the rotor {X}, {Y} and 

the gearbox {X b >, {Y fa } at the bearing locations (Choy 1987). The mass-inertia 
and gyroscopic effects are incorporated in the mass matrix [M] and the gyro- 
scopic matrices [G y ] and [G^] . The coupled torsional equations of motion for 
the single gear-rotor system can be written as: 
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(3) 


In Eq. (3), {F T (t)> represents the externally applied torque and { F Gt ( t )> 

represents the gear mesh induced moment. Note that Eqs . (1) to (3) repeat for 

each single gear-rotor system. The gear mesh forces couple the force equa- 
tions of each gear-rotor system to each other as well as the torsional equa- 
tions to the lateral equations (Choy, 1989; Cornell, 1981; David, 1987 and 
1988) . In addition, there are equations of motion for the gearbox which 
couple the various gear-shaft-systems through the bearing supports. The gear- 
box equations can be written as: 
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and Z-equation 


-!M( r »- Y -l*M( x »-M-° 


MM * M](M *MW ■ 


(6) 


where F bz (t) i- s an excitation function due to external forces in the axial 
direction. Since the bearing is assumed to be uncoupled in the Z-direction, 
Eq. (6) can be solved independently without considering shaft motion. 


The torsional and lateral vibrations of a single individual rotor and the 
dynamic relationships between each gear-rotor-system are coupled through the 
nonlinear interactions in the gear mesh. The gear mesh stiffness varies in a 
periodic nonlinear pattern with each tooth pass engagement period ( August , 
1986; Cornell, 1981; Savage, 1986) and can be represented by a high order 
polynomial (Cornell, 1981; Boyd, 1987), For the coordinate system as shown 
Fig. 1, the following gear mesh coupling equations are established by equating 
forces and moments (Choy, 1989). For the k th gear of the system, summing 
forces in the X-direction results in: 


’ Gxk ” ^ 

i=l,iA 


tki 


-R .8 . 

ci ci 


- R ck 0 ck + 


( X ci - X ck) 


cos CL 


ki 


( Y ci - Y ck) sin a ki lcos a, 


( 7 ) 


ki 


Summing forces in the Y-direction results in: 
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Summing moments in the Z-direction results in: 
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where n is the number of gears in the system. 

Using the modal expansion approach (Choy, 1987, 1988a, 1989, and 1990), 
the motion of the system can be expressed as: 
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where m is the number of modes used to define each motion. 

Using modal expansion and the orthogonality conditions, (Choy 1989, 1990) 

with the bearing forces due to the base motion expressed on the right hand 
side of the equations the modal equations of motion for the gear-rotor system 
(Choy 1989) can be written as: 

X-Z equation 
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The gear mesh forces and moments can also be ^expressed in the modal form, for 


the k th gear with the gear location at the 1 node, as: 
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where k is the gear-shaft-system number, j is the mode number, and 1 is 
the station location of the gear mesh. 

A set of modal equations of motion can also be written for the gearbox 
(Choy 1989, 1990) as: 
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Similarly, the Y-equation can be written as 
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and the Z-equation as 


M(M ♦MN*l* 2 l(M-[M r ( p “ <t >) 


The procedure involves the solution of the coupled modal equations of 
motion between the (gear) rotors and the gearbox structure. The coupling 
effects of the gear mesh and the bearing supports are also expressed in modal 
coordinates such that the global equations are solved simultaneously in modal 
form. A set of initial conditions for both displacement and velocity of the 
global system are calculated from the steady state conditions at the rotor- 
bearing systems and zero vibration at the gearbox. The modal accelerations A, 
B, A t , A fa , B fa , and of the system can be evaluated (Eqs. (11), (13), (15) , 

and (19) to (21)). A variable time stepping integration scheme (the Newmark- 
Beta Method) is used to integrate the modal acceleration to evaluate the modal 
velocities and displacements at the next time step. A regular time interval 
of 200 points per shaft revolution is used in this study except for a refined 
region of smaller steps at the transition from single to multiple tooth con- 
tact. The modal acceleration, velocity, and displacement calculated from the 
transient integration scheme are transformed back into the generalized coor- 
dinates (Eq. (10)). The nonlinear bearing forces and gear mesh forces can be 
evaluated from the velocity and displacement differentials between the rotors 
and the gearbox structure. 


EXPERIMENTAL STUDY 

The gear noise rig (Fig. 2) was used to measure the vibration, dynamic 
load, and noise of a geared transmission. The rig features a simple gearbox 
(Fig. 3) containing a pair of parallel axis gears supported by rolling element 
bearings. A 150-kW (200-hp) variable-speed electric motor powers the rig at 
one end, and an eddy-current dynamometer applies power-absorbing torque at the 
other end. The test gear parameters are shown in Table I. 

Two phases of experiments were performed on the gearbox; (1) static modal 
analysis of the gearbox and (2) dynamic vibration measurements during opera- 
tion. Modal parameters, such as system natural frequenceis and their corre- 
sponding mode shapes, were obtained through transfer function measurements 
using a two channel dynamic signal analyzer and modal analysis software. For 
this experiment, 116 nodal points (with 3° of freedom each) were selected on 
the gearbox housing. Vibration data was collected from accelerometers placed 
on the gearbox at a few of the node points used in the modal survey. Water- 
fall plots of frequency spectra from node 21 on the gearbox top and from node 
40 on the gearbox end are presented in Figs. 4 and 5. The modal frequencies 
(from Figs. 6 and 7) and the gear mesh frequency are shown on figures 4 and 5. 


DISCUSSION 

The experimentally obtained modes, shown in Figs. 6 and 7, represent the 
major vibration modes of the gear noise rig in the 0 to 3 kHz region. 

Although these modal frequencies are only a small part of the total modes of 
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the system, they represent a major part of the total global vibration of the 
system. In order to produce a compatible analytical simulation of the test 
apparatus, a similar set of modes were generated using a finite element model 
of the gearbox structure. Out of a total of 25 modes existing in the analyti- 
cal model in the 0 to 3 kHz frequency region, the eight dominant modes were 
used to represent the gearbox dynamic characteristics. These analytically 
simulated modes are shown in Figs. 8 and 9. As shown in Table II, the natural 
frequencies of the simulated modes are within 5 percent of the measured modes. 
The three-dimensional analytical mode shapes are very similar to the experi- 
mental modes shapes (Figs. 6 and 7). The correlation in the results between 
the analytical model and the experimental measurements confirms the accuracy 
of the dynamic representation of the test gearbox using only a limited amount 
of modes. 

In the experimental vibration spectra of Fig. 4 (from the gearbox top), 
there are prominent peaks at the gear mesh frequency on the traces on 1500 rpm 
(700 Hz), 3750 rpm (1750 Hz), and 5500 rpm (2567 Hz). These frequencies are 
near natural frequencies found in the modal survey (658, 1762, and 2536 Hz, 

respectively). At the highest speed (5750 rpm), the higher modes (2536, 2722, 
and 2962 Hz) are excited by the gear mesh vibration and its sidebands. 

Similar behavior is shown in the vibration spectra from the side of the gear- 
box (Fig. 5 ) . * In both of these figures, the amplitude of the peaks increases 
with speed as the mass imbalance force increases with the square of the rota- 
tional speed. 

The vibration response of the gearbox was simulated by the analytical 
method presented in section II. The gearbox vibration modes were calculated 
by finite element methods. The gear-rotor system response was calculated with 
a discrete rotor model using the matrix transfer technique. The gear-rotor 
system was coupled to the gearbox in modal coordinates to solve for the tran- 
sient vibration response of the global system. The vibration response was 
then transformed into the frequency domain through a fast Fourier transforma- 
tion routine. 

The analytical vibration spectra shown in Figs. 10 and 11 simulate the 
experimental results of Figs. 4 and 5. Both the experimental spectra (Figs. 4 
and 5) and the analytical spectra (Figs. 4 and 5) show excitation at the 
natural frequencies at shaft speeds of 1500, 3500, and 5500 rpm. At the 
higher shaft speed (5500 rpm), several modes are excited by the shaft fre- 
quency and its sidebands. Small differences between the analytical and 
experimental spectra may be due to (1) small errors (<5 percent) in the cal- 
culated modes, (2) nonuniformat ies in the gearbox not present in the modal, 

(3) errors in measurement of rotational speeds and vibration amplitudes, and 

(4) limitations in modeling the bearings as a radial stiffness element only. 


CONCLUSIONS 

Analytical and experimental studies of a single stage gear transmission 
system were performed. Results from the analytical simulations show good 
correlation with experimentally obtained data from the test gearbox. Some 
basic conclusions from this study can be summarized as follows: 

1. A modal synthesis approach can be used to simulate the dynamics of a 
gear transmission system. 
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2. The accurate gearbox model can be developed by correlating the modal 
characteristics from experimental study with those from analytical 
simulations . 

3. The proper choice of modes used in the modal synthesis will reduce 
the number of modes required in the analysis, without sacrificing accuracy in 
the numerical solution. 


REFERENCES 

August, R. and Kasuba, R. , ’’Torsional Vibrations and Dynamic Loads in a Basis 
Planetary Gear System, " Journal of Vibration, Acoustics, Stress, and 
Reliability in Design , Vol. 108, No. 3, July 1986, pp. 348-353. 

Boyd, L.S. and Pike, J.A., "Epicyclic Gear Dynamics,” AIAA Journal, Vol. 27, 
No. 5, May 1989, pp. 603-609. 

Choy, F.K. and Li, W.H. , ’’Frequency Component and Modal Synthesis Analysis of 
Large Rotor-Bearing Systems with Base Motion Induced Excitations," 

Journal of the Franklin Institute , Vol. 323, No. 2, 1987, pp. 145-168. 

Choy, F.K., Padovan, J., and Li, W.H., "Rub in High Performance Turbomachin- 
ery: Modeling; Solution Methodology; and Signature Analysis," Mechancial 

Systems and Signal Processing , Vol. 2, No. 2, Apr. 1988, pp. 113-133. 

Choy, F.K., Townsend, D.P., and Oswald, F.B., "Dynamic Analysis of Multimesh- 
Gear Helicopter Transmissions," NASA TP-2789, 1988. 

Choy, F.K., Tu, Y.K., Savage, M., and Townsend, D.P., "Vibration Signature 
Analysis of Multistage Gear Transmission," 1989 International Power 
Transmission and Gearing Conference, 5th , Vol. 1, ASME, New York, 1989, 
pp. 383-390. (Also, NASA TM-101442 ) . 

Choy, F.K., Tu, Y.K., Zakrajsek, J.J., and Townsend, D.P., "Dynamics of Multi- 
stage Gear Transmission with Effects of Gearbox Vibrations," Proceedings 
of the 1990 CSME Mechanical Engineering Forum, Vol. 2, Canadian Society 
for Mechanical Engineering, 1990, pp. 265-270. (Also, NASA TM-103109) . 

Cornell, R.W. , "Compliance and Stress Sensitivity of Spur Gear Teeth," Journal 
of Mechancial Design , Vol. 103, No. 2, Apr. 1981, pp. 447.459. — — 

David, J.W. and Park, N., "The Vibration Problem in Gear Coupled Rotor 
Systems," 11th ASME Vibrations and Noise Conference, Boston, MA, 

Sept. 29, 1987. 

David, J.W., Mitchell, L.D., and Daws, J.W., "Using Transfer Matrices for 
Parametric System Forced Response, " Journal of Vibration, Acoustics, 
Stress and Reliability in Design , Vol. 109, No. 4, Oct. 1987, ~ ~ 

pp. 356-360. 

Kahraman, A., Ozguven, H.N., Houser, D.R., and Zakrajsek, J.J., "Dynamic 

Analysis of Geared Rotors by Finite Element," NASA TM-102349, AVSCOM- 
TM-89-C-006 , 1990. 


10 


Lim f T.C., Signh, R. , and Zakrajsek, J.J., ’’Modal Analysis of Gear Housing and 
Mounts," 7th International Modal Analysis Conference , Vol. 2, Society for 
Experimental Mechanics, Bethel, CT, 1990, pp. 1072-1078. (Also, NASA 
TM-101445) . 

Mitchell, L.D., and David, J.W. , "Proposed Solution Methodology for the 

Dynamically Coupled Nonlinear Geared Rotor Mechanics Equations," Journal 
of Vibration, Acoustics, Stress, and Reliability in Design , Vol. 107, 

No. 1, Jan. 1985, pp. 112-116. 

Ozguven, H.N. and Houser, D.R., "Mathematical Models Used in Gear Dynamics - A 
Review,” Journal of Sound and Vibration , Vol. 121, No. 3, Mar. 22, 1988, 
pp. 383-411. 


11 



TABLE I. 


TEST GEAR PARAMETERS 


Gear type 

Number teeth 

Module, mm (diametrial pitch in. ) . . . 

Face width, mm (in.) 

Pressure angle, deg 

Theoretical contact ratio 

Driver modification amount, mm (in.) . . 

Driven modification amount, mm (in.) . . 

Driver modification start, deg . . . . . 

Driven modification start, deg 

Tooth-root radius, mm (in.) 

Gear quality 

Nominal (100 percent) torque, N-m(in.-lb) 


Standard involute, full-depth tooth 

28 

3.174(8) 

6.35(0.25) 

20 

1.64 

0.023(0.0009) 

0.025(0.0010) 

24 

24 

1.35(0.053) 

AGMA class 13 

71.77(635.25) 


TABLE II. - COMPARISON OF EXPERIMENTAL 
MEASURED AND ANALYTICAL MODELED NATURAL 


FREQUENCIES 

Experimental , 

Analytical , 

Difference 

Hz 

Hz 

percent 

658 

658 

0 

1049 

1006 

-4.1 

1710 

1762 

3.0 

2000 

2051 

2.6 

2276 

2336 

2.6 

2536 

2536 

0 

2722 

2752 

1.1 

2962 

3012 

1.7 









Dynamometer-, 



Figure 2. — Picture of gear noise rig. 


Node 40- Y 
(on back side) -> 



Figure 3. — Test gear box 
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Mode: 2 

Frequency, 1048.56 Hz 




Mode: 3 

Frequency, 1709.77 Hz 


Mode: 4 

Frequency, 1999.95 Hz 


Figure 6. — Gear box experimental mode shapes. 
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Mode: 5 

Frequency, 2275.69 Hz 





Mode: 8 

Frequency, 2961,71 Hz 


Figure 7. — Goar box oxporimontal modo shapos. 
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Figure 8 — Gear box analytical mode shapes. 
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Mode: 5 

Frequency * 2336 



Mode: 6 

Frequency - 2636 



Mode: 7 

Frequency = 2752 



Mode: 8 

Frequency * 3012 


Figure 9. — Gear box analytical mode shapes. 
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Mesh frequency 



Figure 10 — Analytical vibration frequency spectrum at node 21 (note: The 
prediced modal frequencies are shown at top of figure). 



0 50 / 100 150 200 250 300 350 400 

Mesh frequency-' 


Frequency, kHz 


Figure 1 1 . — Analytical vibration frequency spectrum at node 40 (note: The 
prediced modal frequencies are shown at top of figure). 
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